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We present a general strategy to derive entanglement criteria which consists in performing a 
mapping from qudits to qubits that preserves the separability of the parties and SU(2) rotational 
invariance. Consequently, it is possible to apply the well known positive partial transpose criterion to 
reveal the existence of quantum correlations between qudits. We discuss some examples of entangled 
states that are detected using the proposed strategy. Finally, we demonstrate, using our scheme, how 
some variance based entanglement witnesses can be generalized from qubits to higher dimensional 
spin systems. 
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A necessary and sufficient condition to assert the sep¬ 
arability of a given general quantum state is an open 
problem. For bipartite quantum systems formed by two 
2-level systems (qubits) or a 2-level system and a 3-level 
system (qutrit), the Peres-Horodecki criterion [HU]: that 
is the positivity of the partial transpose (PPT) of the 
quantum system density matrix, provides such a condi¬ 
tion for separability. We still lack such a criterion to 
fully characterize separability in higher dimensional or 
in multipartite quantum systems. Although the PPT 
test constitutes a sufficient condition for detecting bipar¬ 
tite entanglement in the case of systems of arbitrary di¬ 
mension, it requires not only local manipulation of each 
party, but also a full reconstruction of the state density 
matrix. Such requirements may be prohibitive from an 
experimental perspective. This is why entanglement wit¬ 
nesses [ 3 ], which are sufficient criteria for detecting bi¬ 
partite or even multipartite entanglement based on less 
demanding measurements, have attracted so much inter¬ 
est lately [1]. 


In this Letter, we present a general scheme to detect 
entanglement in systems of arbitrary (finite) dimension 
based on the mapping of qudits to qubits. The proposed 
mapping, which actually constitutes a general and op¬ 
erational formulation for dichotomization, preserves the 
separability of the subsystems, ensuring that it does not 
create entanglement that did not exist in the original 
system. Therefore, if the mapped qubits system is en¬ 
tangled, we can assert that the corresponding qudits are 
also entangled. In this way, the proposed mapping en¬ 
ables the application of entanglement criteria originally 
derived for qubit systems to qudit ones. 


We start by defining the mapping Mu oi a, density op¬ 
erator p acting on the d-dimensional Hilbert space to 
a density operator a = Mu(p) acting on a 2-dimensional 


Hilbert space 
Mu : ^ 

p^Mu{p) = Tr ,3 [C/ (p 0 |0) (0|) U^] . 

( 1 ) 

|0) is an arbitrary fiducial state in U is a linear 

operator from the tensor product space to 

another Hilbert space with D > d, and Tr^ 

means a partial trace over the part only. Bi^H) de¬ 
notes the space of the bounded operators acting on "H. To 
preserve hermiticity, positivity and trace of the density 
matrix, it is enough to require that U is an isometry, that 
iaWU = l 2 d , where l 2 d is the identity on [S]. 

Since we are interested in detecting bipartite entangle¬ 
ment in the state described by p acting on 
we define a mapping MuiU^ s-s follows: 

Mu^u2 ■■ ® -7 0 

p ^ r = MuiU^ip): (2) 


with: 


Mut_u-2{p) = TrijjCa 


Hi 0 C/2 (p 0 |0)(0| 0 |0)(0|) U\ 0 Ul 


( 3 ) 

where Ui is an isometric operator acting from 
to(for z = 1,2) andTr^jj D 2 means a partial 
trace over states. Analogously to the case 

with Mu: the mapping MuiU^ preserves positivity and 
trace of the density operator. Furthermore, it also pre¬ 
serves separability: indeed, if p is a 2-qudit product state 
p = pi0p2, thenr = MuiU^ip) = MuApi)®^U 2 {P 2 ) is 
also a product state. This property can be extended to all 
separable states, pure or mixed, by convexity. Therefore 
if r is an entangled state, we are sure that the original 
2-qudit state p G 0 is entangled. As all 

entangled states in 0 have one nonpositive 

eigenvalue for the partially transposed matrix, two nat¬ 
ural questions arise: given an entangled density matrix 




2 


r, what type of entanglement is present in the original 
2-qudit state defined in ® ? And how 

does the set of detectable entangled states depend on 
the isometries Ui and U 2 that are used to implement the 
mapping? 

In the present Letter, we provide some answers to these 
questions in the operationally relevant case where the 
isometries Ui ^2 correspond to a mapping which preserves 
the SU(2) rotational invariance and, in addition, lead 
to an entanglement witness which can be easily imple¬ 
mented experimentally. 

We start by remarking that the isometry U which char¬ 
acterizes the mapping Mu can be parametrized by 4 lin¬ 
ear operators Ai {i = 0,1, 2, 3) from to as fol¬ 
lows: 


shown to be [5]: 

VJTt) = 271 ■ VTTt) ’ 
. 1/6 \ ^ f b b^ \ 

'“272It? v 7 +tJ’ '“272It?^T f+iJ■ 

( 8 ) 

We emphasize that this mapping is practical, as the 3 
expectation values {Ji)p can be easily measured. In ad¬ 
dition, it conserves the rotational invariance. Indeed, 
suppose that we perform a rotation Tifi{a) by an angle 
a around a given vector n. Then, p is transformed as 
7 = ". It is not difficult to show that p' 

is mapped to the rotated qubit : 


3 

U = ^Ai^a^, (4) 

i^O 

where tr/s, for i = 1,2,3, are the Pauli matrices and 
CTo = 1 is the identity operator in The isometric 

property WU = 1 implies that: 


• 3) ^ Mu ■ (9) 

The invariance displayed in Eq. 0 is a consequence of 
the simple vectorial relations: 


{nz\a) \j] )p-a= {J)p, ■a={J)- 7^a(a) [a], (10) 


-A + Aq^Aq = 1, 

Ai Ao -f Ao^A -k iAi A A = 0. (5) 

With the parametrization Eq. 0 the mapping Eq. 0 
can be written as: 


Mu[p)='^{Aj^Ai)pai\Q){Q\aj. ( 6 ) 

i,j=0 

We proceed by introducing the following intuitive and 
convenient specific case for this mapping: 


where J is the vector whose 3 components are the 3 op¬ 


erators 7, for i = 1,2,3, and 


J 


is the corre¬ 


sponding vector in the rotated frame. Since SO(3) and 
SU(2) are isomorphic, we can apply all possible unitaries 
to the mapped qubit by rotating the original angular mo¬ 
mentum correspondingly. 

Now, we consider the case of a 2-qudit state p in 
S(T-L(M (g) and use the mapping Muuj^, where 

each Ui(i = 1,2) implements a mapping as the one given 
by Eq. 0. The resulting mapping MuiU 2 can then be 
written explicitly as follows: 


■^C/(P) = ^ 1^1+j^(7)pcr,j , (7) 

where Ji{i = 1,2,3) are the 3 Cartesian angular mo¬ 
mentum components. They are the generators of SU(2) 
rotations around x, y, and z axis, in the d-dimensional 
Hilbert space The denominator j is the largest 

eigenvalue of 7, i.e., it is such that d = 2j + 1. Einally, 
{Ji)p — Tr[/oJi] denotes the expectation value of Ji. 

The mapping Mu given by Eq. 0 can be shown to 
be a valid mapping defined from Eq. 0, by exhibit¬ 
ing the corresponding A^ operators of Eq. 0. We note 
that the Ai operators can be expressed in a simple form 
with the help of the two bosonic annihilation opera¬ 
tors a and b corresponding to the Schwinger represen¬ 
tation: J+ = aT, J_ = b^a, where J± = Ji ± iJ 2 , and 
J 3 = — b^b) with the restriction aT -I- b^b = 2j. 

By a straightforward calculation the operators Ai realiz¬ 
ing the isometry U through Eqs. 0 and Eqs. 0 can be 


MuiU2{p) = ^ 




2=1 








di=i 


( 11 ) 


where is the i-th (? = 1,2,3) angular momentum 
component on the dfe-dimensional Hilbert space 
{k = 1, 2), with dk = 27 '^'=) -|- 1. 

An important property of this particular mapping is 
that the partial transpose (PT) of a 2-qudit state is 
mapped to the PT of its corresponding 2-qubit state, 
i.e., MuiU 2 iP^^) = A4^7l^72(p)^^ where is the PT 
with respect to the second party. This follows directly 
from the fact Tr[A^^il] = Tr[Ai?^^] for A and B acting 
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on ® and that the transpose of the compo¬ 

nents of the angular momentum operator fulfills = Ji, 
Jj = —J 2 , and Jj = J 3 . As a result, taking into ac¬ 
count positivity preservation, a 2 -qudit state which re¬ 
mains positive after a partial tranpose (a PPT state) is 
mapped to a PPT 2-qubit state [7]. Moreover, it leads 
to an operationally easy to implement entanglement wit¬ 
ness based on second order correlations. This is achieved 
by using the following substitutions: 

{a,^l)Ma,u,(p) = - JIT) -= 




(1® 4^^=)) 


0 (^ 2 ) 


^ rl (12) 


® MuiU2(p) - j(d2)j{d2) - 


which are direct consequences of Eq. ©■ Therefore we 
consider the following form for the general mapped state: 


MuiU^ip) = ^ 


- ( 


■ 7=^.1 ( 




2=1 


(13) 

where we assume that the rotations needed to diagonalize 
Tij have been performed. We have introduced the vectors 
(j = 1 , 2 ), which have as components the (i = 
1,2,3) defined in Eq. (12 1 , after the mentioned needed 
rotations. 

Once we have the form of Eq. 0 or Eq. ([I^, the 


best entanglement witness that we can use is the Peres- 
Horodecki criterion [5] relying on the positiveness of the 
partial transpose. This criterion can be simplified by 
considering the geometric picture which was developed 
in Ref. [8]. In this paper, it was shown that for the 
states of the form given by Eq. (13), the vector T = 
{Til,r 22 ) 233 } G must lie within a tetrahedron with 
vertices {(-1, -1,-1), (-1,1,1), (1, -1,1), (1,1, -1)}, to 
fulfill the positiveness requirement. Each of the 4 
vertices of this tetrahedron is reached when the 2- 
qubit state is one of the four Bell-states 1$^) = 

1212^111, |\ji±^ = jn tJiis picture, the sep¬ 

arable states are those for which the vector T = 
{Til, T 22 , T 33 } G lies within the octahedron with ver¬ 
tices {(±1, 0,0), (0, ±1,0), (0,0, ±1)}. This last property 
can be put in the following more compact form: for any 


separable state of the form given by Eq. (13), T verifies 
the inequality: 


IT11I + IT22I + IT33I < 1 . 


(14) 


Using the definitions given by Eq. (12), Eq. (14) can be 


re-expressed as a 2 -qudit entanglement criterion: 

For any separable state p acting on 0 the 

/p : i = 1 , 2 ,3} verifies the following 


vector 0 


inequality: 




{J 2 

(4 


(A) 


t{^2)\ 

^2 Ip 




Therefore, all states that violate inequality (15) lie out¬ 


side the octahedron and are thus entangled. This 2-qudit 
entanglement criterion has the advantage of being very 
simple to test experimentally. One can notice that there 
are no 2 -qudit states that are mapped to any of the 2 - 
qubit Bell states [5]. Therefore, the vertices of the tetra¬ 
hedron do not belong to the image of our mapping. 

To have an insight about the efficiency of our crite¬ 
rion to detect entanglement, we apply it to two known 
families of qudit states that have been extensively stud¬ 
ied in Refs. [HHia. From now on, we suppose that the 
two qudits are of the same dimension, that is di = d 2 = 
d = 2j + 1. First, we recall the family of d^ maximally 
entangled 2 -qudit states |Ufc;)(A:, ^ = 0 , 1 , • • • , d — 1 ) that 
generalize the four 2-qubit Bell states [SHU]: 


1 


d-1 


|f2/c/) = Wki 0 l|Uoo), with |Uoo) = ^ ^ |m,TO),(16) 

where the W^i operators acting on the first qudit are 
the Weyl operators defined as 

Wki\m) = e a ^ |(to - Z)mod d)- (17) 

These 2-qudit Bell states Pki = |Ufc;)(Ufe/| are locally 
maximally mixed states, that is, by taking their partial 
trace one obtains the maximally mixed state A It is 


interesting to notice that they are mapped by Eq. (13) 


to a locally maximally mixed 2-qubit state. Indeed, we 
have: 


j 


and 


"kl 


(1 0 j) 


Pki 


= 0 . 


For such states, our simple criterion Eq. (15) is as strong 


as the PPT criterion applied to states given by Eq. (13) 


and detects all entangled locally maximally mixed states 
that are detected by the latter. We can thus say that ev¬ 
ery maximally entangled 2 -qudit state \flki) is detected 


by our criterion Eq. (15), even though these states are 
not mapped to 2-qubit Bell states [S]. Instead, they 
are mapped to mixed states that are locally maximally 
mixed, that is, a convex sum (statistical mixture) of 2 - 
qubit Bell states. 

We proceed by exploring the statistical mixtures of 
maximally entangled 2 -qudit pure states that can be de¬ 
tected by our criterion Eq. (15). We start by applying 


our criterion to the so called Werner states which form a 
good description of the effects of phase and depolarizing 
noise in maximally entangled states nniisi: 


Pa — Q!|Uoo)(Uoo|- 


1-0 -1 

_i, „th _ 


< o < 1 . (18) 
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The bounds for the parameter a are such that pa is posi¬ 
tive. It is known [13] that pa is entangled iff < a < 1. 
A straightforward application of our criterion Eq. (15) 
to Pa gives that if aj (j -I- 1) > then pa is entan¬ 
gled. Therefore, our criterion can detect all the entangled 
states Pa for a G [ {j+i) i !]■ Recalling that d = 2j -|- 1, we 

realize that entangled states with 2 {l+i) < ® are 

not detected. 

As a more specific example, we now consider the 3- 
parameter family of 2 -qudit states: 

_ l-a-/3-7 

Pcx,0,-y — /n / I -I \9 1 H” Q^-^00 




(2j + 1)2 

2 i 

'^P^0 + 


7 

2 j + l 


2 i 

2=0 


(19) 


where the Pki are the projectors on the \ flki) states. This 
family of states is a generalization to arbitrary dimen¬ 
sional qudits of the 2 -qutrit states originally introduced 
in Refs. dOHIH to study bound entanglement. Density 
matrices as given in Eq. (371 are interesting because their 


eigenvalues and those of their partial transpose can be 
explicitly expressed as a function of the parameters a, 13 
and 7 | 6 |. This allows to locate the set of PPT Pa,i 3 ,-y 
in the space spanned by a, /3, and 7 . To ensure positiv¬ 
ity, the parameters a, /3, and 7 , must verify the following 
inequalities |H]: 

l-a-/3-7^^ , 1-Q!-/3-7^^ 

-“ •“+ (2J + 1F 


1 + 

(2j 


{2j + 1)2 

1 — a — /I — 7 


1)2 


>0 + 

- ’ {23 + 1? -"• 


These inequalities define the interior of a tetrahedron. 


Now, applying our criterion Eq. (15), we obtain that all 
separable states Pa,p,'y given by Eq. (19) are such that: 


\{j + l)(a + I3) + {j - 2 ) 7 ! -t {j + 1) |2a - /3/j| 
3j 


< 1 . 


. 

All states Pa,p,-y for which the inequality Eq (38) is vio¬ 
lated are entangled. 

In order to have an idea about the efficiency of crite¬ 


rion Eq. (15), we have calculated the eigenvalues of the 
PT of Pa,/ 3,7 and thus obtained explicit conditions on 
q:,/ 3, and 7 for Pa^^^ to be positive (Sj. In Fig. we 
present the tetrahedron of positive states Pa,i 3 ,'y given by 
Eq. (37) in parameter space for the cases j = 1 (Fig. [3 a) 
and j = 2 (Fig. 0 b)- The set of PPT states are depicted 
in blue, so the remaining volume of the tetrahedron corre¬ 
sponds to entangled states. We have represented in red 


the states that are detected by our criterion Eq. (15). 


We see clearly that it detects a significative part of the 


entangled states parametrized by Eq. (37). Nevertheless, 


comparing Fig.[2a and Fig.[^b, we note that this volume 
decreases when j increases from j = 1 to j = 2. However, 



FIG. 1. Geometrical representation of states Pa,/ 3 ,.y [Eq. (37l] 
in parameter space for (a) -|- 1 = 3 and (b) -|- 1 = 5. All 

physical states lie inside the tetrahedron whereas the PPT 
states lie inside the blue region [a cone for the 2-qutrit case 
(a)]. Red regions depict entangled states detected using our 
criterion Eq. (151, whereas the yellow region hosts the non 
detected entangled states by this criterion. 


different criteria which present a better scaling with di¬ 
mension for this particular family of states can be found 
by changing the isometry U in Eq. @ (or equivalently 
changing the corresponding Ai operators in Eq. , used 
to map each qudit to a qubit). 

Until now, we only have considered the entanglement 
of two qudits. Now, we adress the problem of detecting 
entanglement in a large N qudits system. An interest¬ 
ing consequence of our mapping is that it can be easily 
extended to map a system of N qudits to a system of N 
qubits. Indeed, by applying Eq. ([^ individually to each 
qudit, the separability among the parties is preserved. If 
we denote La = ^{a = 1, 2, 3) and Lq = 1, then the N- 
qubit mapped state corresponding to the iV-qudit state 
p can be written as 


^®^u{p) — E {®f=lLki) p®f=i(^ki-{2P) 

fce{0,l,2,3}« 


An immediate consequence is that for any k G 
{ 0 , 1 , 2 ,3}^, we have 


(®i=l^fci)p — {®i=l^ki)M^N„{p)- (22) 


By using this property to compute first and second order 
correlations, we can provide an alternative derivation of 
the spin squeezing inequalities detecting Wqudit entan¬ 
glement from the fV-qubit ones introduced in Ref. [T3] . It 
was shown in Ref. m that all separable iV-qubit states 
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satisfy the following inequalities: 


{AAf + {AS2f + iAS3f>^, 


( 5 ^) + (S^P) -(N- l)iApr < 
(5^) -{N- l)[{ASpf + (A4)2] < 


(23) 


where Sa = ^ ®’a is the collective spin operator 

in direction a. The indexes (a, (3, 7 ) may assume any 
permutation of (1,2,3) and the following definitions have 
been used: 

(A5„)2 ^ (52) - (5„)2, 

(5^) = (24) 

(A5„)2 = {Sp)-{S^)\ 


Using Eq. (22 1 , we obtain, starting from the A-qubit in¬ 
equalities Eqs. (23), the inequalities satisfied by all sep¬ 
arable A-qudit states (A spins j). This is achieved with 
the following substitutions HU: 


5„ ^ ^ Ja, 5^ ^ (25) 

where and = q E^i=i 4 ® P- 

Therefore, the relation between the entanglement crite¬ 
rion for A-qubit and A-qudit systems, which was already 
considered in Ref. mi, can be thought as a simple con¬ 
sequence of the particular mapping explicitly given by 
Eqs. ^ or equivalently by Eq. 0- As a consequence, 
we have shown that the qudit entanglement revealed by 
the qudit spin squeezing inequalities can always be recast 
as qubit-like, or dichotomic, spin squeezing. Thus, qudit 
spin squeezing inequalities do not evidence entanglement 
of higher dimension than the qubit squeezing ones. Fi¬ 
nally, we notice that by choosing different A, operators 
in Eq. (|^ we can expect to find new multipartite qudits 
entanglement criteria. 

In conclusion, we have presented a general scheme to 
map qudits to qubits that can be used to define crite¬ 
ria to detect some type of entanglement between qudits. 
We have applied this general scheme to provide a spe¬ 
cific entanglement criterion based on the measurement 
of qudit-qudit correlations. 

In addition, our results provide a way to classify multi¬ 
partite qudit entanglement according to its detectability 
through dichotomization. Finally, it opens the interest¬ 
ing question of what are the specific entanglement types, 
if any other than bound entanglement, that fail to be 
detected by our method. 
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SUPPLEMENTAL MATERIAL 

This supplementary material contains 3 parts: In the 
first part, we prove Eq. ( 8 ) of the main text. In the sec¬ 
ond part we prove that the 2-qubit Bell states are not in 
the image of our 2-qudit mapping given by Eq. (11). Fi¬ 
nally, in the third part we introduce the three-parameter 
family of 2 -qudit states that generalizes the family of 2 - 
qutrit states defined in Ref. [ 8 ]. We show how we obtain 
the separability criterion Eq. (20) of the main text as 
a function of the three parameters defining these states. 
Next, we calculate the eigenvalues of these states and 
their partial transpose. 


PROOF OF EQ. (8) 

In accordance with our general scheme for mapping a 
qudit state to a qubit one defined in the text, we need to 
find operators Ai defining the isometry U = 
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such that 


ip) — 'y ' i^j^^i)p 

jj =0 


(26) 


is equal to Eq. (7) of main text which we recall here: 

3 \ 


Muip) = \[t + -.Y.^J^) 


1 


3 


(27) 


In order for U to be an isometry WU = 1, we have the 
following conditions: 


^ 0^0 “ 1 “ 71 | Ai + ^ 2^12 + ^ 3^3 “ 1 : 

“t” ^ 0^1 ^^2^3 — ^^3^2 — 0 i 

“t” -^0^2 — ^^1^3 ~t“ iA^Ai = 0, 

+ -Agj43 + ^^1^2 — iA^Ai = 0, 


(28) 


Eq. (11) of the main text. This can be proved by con¬ 
tradiction: consider that there exists a state p such that 
AiuiU^ip) is one of the Bell states, say |4'+)(^'+|. For 
this state we have: 

Tr[|^'+)(^+|a3] = 1 , 

Tr[|T+)(^+|a 2 ] =-1, Tr [|4-+)(4-+|ai] =1. (30) 

From Eq. ( [M| and Eq. (11) of the main text, we find 

-vvtiich in turn implies that the 
state p must be a pure state of the form + 

f 3 \—j(d-i)^—jid 2 )'j with \af + |/3|^ = 1. For such states 
the other values of {o'i)MuiU 2 ip) for z = 1 or * = 2 are 

= (J^) 0 


'IU 1 U 2 (p) 

zero and not 1 or - 1 , indeed {j[‘^ 

M2)^P 

can be made for each Bell state. 


JT^’)p = 0 for di > 2 or ^2 > 2. The same reasoning 


13^U ^ 0 ^ 

wich are equivalent to Eq. (5) of main text. Comparing 
Eq. (|27|) and Eq. ([2^, we get 


— 2242.^3 “t“ ^-^ 3.^2 T 

2241^43 — iA^Ai -\- 


^1^3 “t“ 242^41 2^gA2 — 2 j4224q 


J 2 I 


3-PARAMETER FAMILY OF 2-QUDIT STATES 
Preliminaries 


^2^3 4“ 7 I 2 A 2 — 2 AqAi 4“ 2 A|Ao *^2 [' 29 ') First, we recall the family of cP maximally entangled 

2 + 2 ^ 


1 J 3 


— 2AJA2 4- 2^2^! 4- ^ 3^3 + ^0^0 ~ 2 ~ 2j 

If we define ^03 = v7 (^0 4-A 3 ) andAi 2 = v(J (Ai 4 - 2 ^ 2 ), 
the set of equations Eqs. (29) can be simplified using the 
set of conditions Eqs. @10 


2 -qudit states |Ofc/)(fc, Z = 0 , 1 , • • • d — 1 ) that generalizes 
the four 2-qubit Bell states as introduced in Ref. m- 


AJ 3 A 03 — 


3 _ J-’^3 A _ 1 + ^-^2 

, ..a1j^ 2“^12 — „ j -^03“^12 “ o 


These equations can be simplified with the help of the 
two bosonic annihilation operators a and b corresponding 
to the Schwinger representation of the qudit; J+ = a^b, 
J- = b^a and J 3 = ^( 0^0 — b^b) with the restriction 
a'^a + b'^b = 2j where J± = Ji ±iJ 2 - We finally get: 


d-l 

\^ki) = Wki 0 l|lloo) with |Ooo) = ^ E I^>"^)(i31) 

where the cP Wki operators acting on the first qudit are 
the Weyl operators defined as: 

Wki\m) = 22 ;'"(”"“'^|(m - Z)inod d), where w = e^’"*/‘^(32) 

These operators verify the following orthogonality rela¬ 
tion in the Hilbert-Schmidt norm: 


Tr 


AJ3A03 — o J AI2A12 — „ , AQ3A12 — „ . 


T 


T 


WliW^n 


— d Skm^ln- 


(33) 


It is straightforward to show that the operators Ai de¬ 
fined by Eq. ( 8 ) of the main text verify the conditions 
of Eqs. 


and Eqs. (29). Thus, we have proved the 
validity of mapping ([^|. We note that the requirement 
for U to be an isometry and not a unitary operator is 
crucial in the present case. Indeed, the operators Ai do 
not fulfill the conditions for U to be unitary. 


BELL STATES ARE NOT REACHED 

We proceed to show that none of the four 2-qubit Bell 
states belongs to the image of the mapping given by 


These 2-qudit Bell states Pki = |Ofc;)(nfe/| are locally 
maximally mixed state, that is, their partial trace gives 
the maximally mixed state |. 

Using the following relation m- 


d-l 




(34) 


1=0 


each state Pki can be written in the basis Wki 0 W„ 
follows miin]: 


as 


Pki = 


dP 


d-l 

E 

m,n=0 


^ml-knwl^_^ 0 Wi_„. 


(35) 
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The operators J 3 , J+, J_ can also be written in the Weyl 
basis: 


2 j 

1=0 

2i 2j 

j+^Yl Vpw~‘wi,-i, j- = Y^ ipWi^i, 

1^0 1^0 


(36) 


where we have defined 

2 j 


m—0 
2 j 

£ 




m =0 


2 j + l 

\/j(j + 1 ) - {m- j) {m- j + 1 ) _ 

2j + 1 


ml 


From Eqs. (35), (36), and (33), we get 


2 i 


(J3 ® ■^ 3 )|nfc,) = Ti' [PkiJa -/a] = ^ 




m—0 


2i 


(J+ (g) J+)|nfci) = Tr [PkiJ+ g J+] = w ^ ?7™?7i 


m -m ml 
P ^ 5 


m =0 


(J+ g = Tr [PkiJ+ g J-] = 0. 


Now, we have the tools that enable us to study the three- 
parameter family of 2 -qudit states with d = 2 j-|-l defined 
in the main text: 

_ l-a-/3-7 

da,/3,7 — /n _■ 1 -1 \9 H “t“ O^foQ 


(2j + 1)2 
2 i 
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2 i 




2j 21 ^ 2 j -f 1 21 - 

i=l •’ i=0 

For these states we have 

jU + 1) (a + P) + j{j - 2)j 


(37) 


(J 3 ® >^ 3 ) 

(Jl g Jl)pa,p,y - 
(J 2 g = 

so that the criterion 


j{j -fl) (q - P/2j) 

3 

j{j -f 1) (a - /3/2j) 




+ 




(dl) Ad2) 


< 


introduced in the main text for separable states Pa,/ 3 ,~f 
can be brought to the following form: 


|(j -f l)(a + P) + {j - 2 ) 7 ! -f (j -f 1) |2q; - P/j\ 

3j 


< 1- (38) 


Positivity and partial transpose of pa,p,~, 


From the definition in Eq. (37), Pa,p,-y is already in its 
diagonal form and we can easily identify its eigenvalues 

1 — (X — 13 — 7 

and eigenvectors. The eigenvalues are:— 7 —- —15 - 1 - 

( 2 j + 1 )"^ 

^ with degeneracy 2 j -|- 1 , — ^ ^ ^ + — with 


2 j + l 


( 2 j + 1)2 2 j 


1 — (X — [3 — 7 

degeneracy 2j, ■— 7 —-——■ -I- a with degeneracy 1, 

\2j + 1 ) 

and ^ ( 2 £i )2 degeneracy (2j -f l)(2j - 1). To 

ensure the positivity of Pa,p,-y, each of these eigenvalues 
must be positive, hence we obtain the positivity condi¬ 
tions stated in the main article: 


1-q:-/3-7 1-q;-/3-7 

(2j + l)2 

> 1 -^-^ + 2j7 ^^^ 


{2j + 1)2 

P 1 — a — /3 — 7 

^ ^ ( 2 j + 1)2 


( 2 j + 1)2 


As for the partial transpose (PT) of Pa,p,-yi it was shown 
in Ref. [T7] that for a state of the form p = Y4ki=o P'kiPki, 
its partial transpose can be written as: 


=0 -.81= 


(39) 


where, the matrix elements of are defined as m- 


1 


2 j 


{s\Bm\t) = *) : s,t = 0, • • • ,2j. 

^ k^O 

(40) 

It was also shown m that for integer values of j, all 
matrices B^ are unitarily equivalent, while for j half¬ 
integer, there are two classes of unitarily equivalent ma¬ 
trices, the class of B^ for even m and the class for odd 

m. 

In what follows, we will restrict ourselves to the case of 
J integer and we will calculate the elements of the matrix 
Bq. The case of half-integer j can be easily obtained 
based on the calculation of the matrix Bq. 


From Eqs. (37) and (40), we find that the matrix Bq 


has the following form: 


fboo 0 • • • 

0 611 0 


Bq = 


0 

K 


0 


0 0 .’ b2j-i,2j-i 0 

\ 0 K 0 • • • 0 b2j,2jj 


(41) 




































with 




= (0|So|0) = 


(i + l|-So|j + 1) — 


1 — a — /? — 7 a + j3 
( 2 j + l )2 + 2 j + l’ 

1 — a — /3 — 7 


( 2 j + l )2 2 j + l’ 


bk,k = (fc|-Bo|fc) = f ^ : fc e {0, •• • ,2j} \ {0, j + 1}, 


K = (s|-Bo|t) = 

(s|Bo|t)=0 


(2j + 1)2 

a - /3/2j 
2 j + l 


: s + t = 0 and s ^ t, 
: s ands + t 7 ^ 0 . 


The characteristic polynomial of Bq can be easily com¬ 
puted which allows to calculate the eigenvalues of Bq. As 
7 i® direct sum of matrices Bm that are unitar- 
ily equivalent to Bq, then the eigenvalues of Bq are also 
those of ^ . We distinguish 2 cases: for j = I and 
j > 1. For j = 1, the eigenvalues are: 

1 -b 2j{a + f3) - 7 

(1 + 2 j )2 ’ 

1 / 7 _ 2{-l + a + P + -f) _ j {-2ja + 13)"^ + 

2\l + 2j (l + 2 j )2 Y j 2 (i + 2 j )2 J’ 

1 / 7 _ 2{-l + a + P + j) j {-2ja + PY + 

2 l^l + 2 j (l + 2 j )2 Y j 2 (i + 2 j )2 j 

with degeneracy 3 each. Whereas for j > 1, the eigen¬ 
values are: 

1 -b 2j(a -b/3) - 7 
(1 + 2jy ’ 

-P - 2j(-l - 2ja -b 2/3 -b 7 ) /3 - 2j(-l + 2(1 -b j)a -b 7 ) 
2 j (1 + 2 jT ’ 2 j(l + 2 j )2 

1 / 7 _ 2(-l-bQ:-b/3 + 7) _ l {-2ja + PY 

2 l^l + 2 j (l + 2 j )2 Y j 2 (l + 2 j )2 J’ 

1 / 7 _ 2(-l-bQ:-b/3 + 7 ) j {-2ja + I3Y 

2 l^l + 2 j (l + 2 j )2 ^Y j 2 (l + 2 j )2 J 

with degeneracy 2j + 1 , ( 2 j -b l)(j - 1 ), ( 2 j -b l)(j - 1 ), 

(2j -b 1), and (2j -b 1) correspondingly. Thus by imposing 
positivity on the above eigenvalues, we get the set of 
conditions for the state Pa,/ 3,7 to be PPT. This is how we 
compute the blue region on figure 1 (a) and (b) in the 
main text. 



